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ABSTRACT 

The Stripe 82 Massive Galaxy Catalog (s82-MGC) is the largest-volume stellar mass-limited sample 
of galaxies beyond 2 ; Ri 0.1 constructed to date. Spanning 139.4 deg^, the S82-MGC includes a mass- 
limited sample of 41,770 galaxies with logM,/M 0 > 11.2 to z « 0.7, sampling a volume of 0.3 Gpc*^, 
roughly equivalent to the volume of the Sloan Digital Sky Survey-I/II (SDSS-I/II) z < 0.15 main 
sample. The catalog is built on three pillars of survey data: the SDSS Stripe 82 Coadd photometry 
which reaches r-band magnitudes of ^23.5 AB, YJHK photometry at depths of 20*^*^ magnitude (AB) 
from the UK Infrared Deep Sky Survey Large Area Survey, and over 70,000 spectroscopic galaxy 
redshifts from the SDSS-I/II and the Baryon Oscillation Spectroscopic Survey. We describe the 
catalog construction and verification, the production of 9-band matched aperture photometry, tests of 
existing and newly estimated photometric redshifts required to supplement spectroscopic redshifts for 
55% of the logM*/M 0 > 11.2 sample, and geometric masking. We provide near-IR based stellar mass 
estimates and compare these to previous estimates. All catalog products are made publicly available. 

The S82-MGC not only addresses previous statistical limitations in high-mass galaxy evolution studies 
but begins tackling inherent data challenges in the coming era of wide-field imaging surveys. 

Subject headings: catalogs — galaxies: abundances — galaxies: general 


1. INTRODUCTION 

A new era of truly panoramic imaging surveys has be¬ 
gun that promises new insights into fundamental ques¬ 
tions in cosmology and galaxy evolution from z < 1 
that have been hindered by the relatively small vol¬ 
umes and resulting statistical uncertainty available to 
date. Adva ncing the legacy o f the Sloan Digital Sky Sur¬ 
vey fSDSS. lYork et al]l200nfl . surveys like DE^, HSCFI. 
KiD ^^. VIKING^3- DECaL^^i and eventually Euclid, 
LSSTTd. and WFIRST will provide thousands if not tens 
of thousands of square degrees of deep multiband imag¬ 
ing data. 

These data sets offer exciting opportunities. Beyond 
cosmological constraints, from the galaxy evolution per¬ 
spective, they will enable precise measurements of evolv- 
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ing number densities that can chart galaxy growth and 
the rate of flow between transitioning populations, thus 
constraining the physical drivers of evolution. They also 
present new challenges, from questions of how to self- 
consistently process and analyze huge data volumes, to 
the rising importance of subtle systematic errors in the 
estimators we use to derive (photometric) redshifts, to¬ 
tal fluxes, intrinsic colors, and physical properties such 
as star formation rate and stellar mass, M*. 

Ahead of the rising tide of “Big Data” in astronomy, 
we can already make progress on “large-volume” ques¬ 
tions and begin tackling some of the challenges above 
using extant surveys. The Canada-France-Hawaii Tele¬ 
scope Legacy Wide Survey (CFHTLS-Wide) is a pioneer¬ 
ing example in the ~100 deg^ regime. The “Stripe 82” 
region, spanning a 2/5 wide by 110 degree long stretch of 
the celestial equator in the Southern Galactic Cap, is an¬ 
other example. As we describe below, Stripe 82 is not as 
deep as CFHTLS-Wide, but offers two features that we 
exploit in this paper to construct the Stripe 82 Massive 
Galaxy Catalog (s82-MGC), the largest-volume near-IR 
selected M»-limited sample of galaxies beyond z « 0.1 
assembled to date. 

First, Stripe 82 was the subject of repeated imag¬ 
ing in SDSS and therefore reaches ugriz depths that 
are roughly 2 magnitudes de eper (90% comple teness 
for galaxies at r ^ 22.5 AB, lAnnis et al.l I2014H than 
the single-epoch imaging. An “SDSS Coadd” of these 
data was processed and analyzed to produce a pub- 
licly available catalog as part of SD SSD at a Release 7 
( Abazaiian et al.l 120091 : lAnnis et al.l 1201411 . iJianer et "m 
(j2014li provide an alternate set of Coadd images, but 
no photometric catalog. The added depth in the Coadd 
is critical for obtaining reliable photometric redshifts 
(photo-zs) for massive galaxies (logM,/M 0 > 11 ) that 
can be used to supplement the color-selected spectro- 
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scopic samples out to z ~ 0.7. Second, enabling more ro¬ 
bust stellar mass estimates and better measures of spec¬ 
tral energy distributions, near-IR photometry in Stripe 
82 is avail able from the UKIRT I nfrared Deep Sky Survey 
(UKIDSS, [Lawrence et al.l[2007ll . Specifically we use the 
YJHK photometry from the UKIDSS Large Area Sur¬ 
vey (LAS) component which reaches a depth of AB^20 
over roughly 230 deg^ in Stripe 82. 

Stripe 82 has also featured prominently in spectro¬ 
scopic campaigns. The SDSS-I/II m ain sample, the 
lumin ous red galaxy (LRG) sample (|Eisenstein et al.l 
1200111 . and the Baryon Osc illation Spectroscopic Survey 
(BOSS lPawson et al.||2013 D as well as 2 SLAQ, 2dF, 6 dF, 
and WiggleZ (|Drinkwater et aP 1201011 all cover Stripe 
82. In addition, the stripe has been observed by nar- 
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1201111 . This array of spectroscopic redshifts is a criti¬ 
cal scientific asset to the S82-MGC. The large number of 
spectroscopic redshifts (spec-zs) from SDSS (especially 
BOSS) provides the foundation needed to build a com¬ 
plete, M*-limited sample. Roughly 45% of galaxies in 
the S82-MGC sample with logM*/M 0 > 11.2 have spec- 
zs. Deeper spec-zs provide additional checks and train¬ 
ing sets for photometric redshifts (photo-zs) that are re¬ 
quired to supplement incomplete spec-z samples in the 
S82-MGC. 

A predecessor to the S82-MGC was present ed in a pio¬ 
neering effort bv iMatsuoka fc Kawaral (120101 1 who down¬ 
loaded and reprocessed the UKIDSS-LAS (DR3) Lf-band 
imaging data and a large fraction of the SDSS Coadd 
ugriz imaging in a 55 deg^ region of overlap in Stripe 
82. The UKIDSS sources were smoothed to a simi¬ 
lar spatial resolution as the S DSS and photometry per ¬ 
formed with Source Extract or (|Bertin fc Arnout^llOOOH . 
iMatsiioka fc Kawaral (|20ir)l l use their catalog to con¬ 
strain the growth history of massive galaxies, finding evi¬ 
dence for a more rapid increase since z 0.9 in the num¬ 
ber of the most massive galaxies {log M^/Mq > 11.5) 
compared to a less massive sample ( 11.0 < logM,/MQ < 
11.5). The construction of the S82-MGC was motivated 
by similar goals which will be the subject of a future pa¬ 
per. In addition to the nearly three times larger area cov¬ 
ered by the S82-MGC and the large number of spec-zs now 
av ailable, a major difference o f our approach compared 
to iMatsuoka fc Kawaral (|2010f) is the use of the synmag 
synthetic aperture photometric matching technique that 
works at the catalog level without requir ing reprocessing 
of imaging data fsee lBundv et al.ll2012l l. With the nec¬ 
essary catalog information in hand, SYNMAG results can 
be obtained quickly for data volumes of nearly any size, 
making 9-band matched photometry over the full Stripe 
82 relatively easy to achieve. 

As the S82-MGC was being assembled, new surveys 
in Stripe 82 started on several facilities: CS82 (the 
CFHT Stripe 82 Survey, Kneib et ah, in preparation), 
an z-band imaging survey to zab ~ 23.5 covering 173 
deg^ with a median seeing of 0'.'6; the VISTA-CFHT 
Stripe 82 Survey (VICS82, Geach et ah, in prepara¬ 
tion), which uses the Visible Infrared Survey Telescope 
for Astronomy (VISTA) and the GFHT Wide-held In¬ 


fraRed Gamera to image 140 deg^ in the J and Ks 
bands to 22 AB (5cr point-source); and the Spitzer- 
IRAG Equatorial Survey (SpIES, PI: Richards), which 
has now mapped 100 deg^ in the Spitzer 3.6^m and 
4.5/iTO channels to 5(t depths of 22.8 AB and 22.1 AB, re¬ 
spectively. Adding information from these observations 
and others to S82-MGC will be valuable. Another valu¬ 
able dataset is the near-IR ph otometry of WISE (Wide- 
held Infrared Survey Explorer I Wright et al.l[ 2 (TT(il l. New 
catalogs o f WISE-bSDSS in atched photometry are now 
available (iLang et al.l 1201411 . The WISE W1 and W2 
channel depths are roughly comparable to the UKIDSS 
LAS. Here we focus on the UKIDSS photometry, espe¬ 
cially for the additional near-IR measurements that are 
useful for constraining galaxy SEDs. 

In this paper, the first in a series, we describe the 
construction of the Stripe 82 Massive Galaxy Gata- 
log, its underlying datasets, various tests and valida¬ 
tion efforts we have undertaken, and value-added prod¬ 
ucts we have derived. The S82-MGC is publicly avail¬ 
able from massivegalaxies.com. The SYNMAG tool 
necessary for fast, catalog-level matched aperture pho¬ 
tomet ry in the S82-MGC was presented in iBundv et al.l 
(|20I2f) an d makes use of Ga ussian Mixtures as de- 
scribed in iHogg fc Land (j2013[ l. Paper II in the series 
(jLeaiithaud et al.l 120151 1 uses the S82-MGC to charac¬ 
terize the M* completeness of the BOSS spectroscopic 
samples. Paper III (Bundy et al. 2016) analyzes mass- 
complete samples drawn from the S82-MGC to set con¬ 
straints on evolution in the galaxy stellar mass func¬ 
tion and growth among the most ma ssive galaxies since 
z ~ 0.7. The S82-MGC is also used in iSaito et al.1 (1201, 5f l 
in their characterization of the host dark matter halos of 
BOSS galaxies and in Ghabonnier et al. (in preparation) 
in a study of highly compact massive galaxies. 

The paper is structured as follows. Section [2] summa¬ 
rizes the SDSS Goadd and UKIDSS imaging as well as 
the spectroscopic samples that form the basis of the S82- 
MGC. The application of SYNMAGs to derive matched 
photometry is described in Section O as is our method¬ 
ology for building new total flux measurements in the 
near-IR given the significant problems in the UKIDSS- 
LAS estimators. Star-galaxy separation is discussed in 
Section m where near-IR colors are used to provide a sig¬ 
nificant improvement over a classification based on SDSS 
shapes alone. We compile and test photometric redshifts 
in Section [51 In Section [5] we analyze the geometry of the 
s82-MGC’s on-sky footprint using the MANGLE software 
package and apply rejection masks. A clean, near-IR 
limited subsample, ukwide, is presented and character¬ 
ized in Section III Section [ 8 ] describes our derived mea¬ 
surements of galaxy properties, including M, estimates, 
which we compare against previous values produced by 
the BOSS team. 

Throu ghout this paper, we use the AB magnitude 
system (lOke fc Gunnlll98'^ and adopt a standard cos¬ 
mology with Hq=70 /i 7 o km s“^ Mpc“^, Dm=0.3 and 
Da=0.7. Stellar mass e stimates assum e a Ghabrier ini¬ 
tial mass function fIMF. lcThabrieill200,3l l unless otherwise 
noted. 

2. SURVEY DATA 

The S82-MGC is built on three pillars of survey data— 
the SDSS Stripe 82 Goadd catalog, the UKIDSS LAS, 
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and the BOSS spectroscopic redshift samples. We sum¬ 
marize each of these in this section. With the SDSS 
Coadd as the starting point, the S82-MGC consists of 
15,342,585 objects, including both extended and point 
sources. After matching to UKIDSS photometry (Sec- 
tion l3.1|l . applying cuts in limiting depths fSection r2.2.1ll . 
improved star galaxy separation (Section|j), and geomet¬ 
ric masking (Section [G]), we derive an M*-limited sub¬ 
sample of S82-MGC that we refer to as ukwide (Sec¬ 
tion [71). The UKWIDE sample consists of 517,714 galax¬ 
ies, with 41,770 above the nominal completeness limit 
oflogM»/M 0 > 11.2. The publicly available (go to 
massivegalaxies. com) input and derived data products 
associated with the S82-MGC are summarized in the Ap¬ 
pendix. 

2 .1. SDSS Stripe 82 Coadd catalog 

The “SDSS Coadd” refers to the stacked imaging and 
photometric cata log in Stripe 82 ( - 50° < q:j 2 ooo < +60°) 
first prese nted inlAbazaiian et ^ (1200911 and further de¬ 
scribed in lAnnis et al.l (|20l3b The Coadd was made 
by processing ^^90 repeated visits to the stripe dur¬ 
ing the SDSS-1/11 imaging campaign, including obser- 
vations obtained du ring the SDSS Supernova Survey 
(|Frieman et al.l 1200811 . The source detection, flagging, 
profile fitting, and ph otometry was perfo rmed on data 
taken before 2005 by JAjmis_et_^Jj fl201^ by applying 
the Photo pipeline (iLunton et al.l 1200^ . The SDSS 
camera is described in iCunn et al.l (|l998f l. the telescope 
system in iCunn et al.l ^20061 1. arid the filter system in 
iFukugita et al.l ( ~ 99611 . The 50% point-source complete¬ 
ness limit of the SDSS Coadd is r ^ 24.4 (AB) and 
the 90% galaxy completeness limit is r ^ 22.5 (AB). 
We used the SDSS Catalog Archive Servei0 and the 
PhotoPrimary view to query all sources in the Coadd 
photometric catalog, Stripe82, with TYPE values of 
3 or 6. The resulting catalog has 15,342,585 sources 
and serves as the basis of the S82-MGC. Please see 
lAnnis et al.l (|2014ll and the SDSS website for further de¬ 
tails on executing the query. 

A number of checks on the C oadd photometry were 
performed by lAnnis et al.l (|2014ll . In building the S82- 
MGC, we also compared the colors and the magnitudes as 
measured in the SDSS Coadd to those from the Canada- 
France-Hawaii Teles cope Tensing Survey (C FHTLenS) 
photometry catalog (iHildebrandt et al]l2012ll for sources 
overlapping the CFHT W4 field. We found excellent 
agreement and no evidence for biased photometry in 
the Coadd. Further support comes from the quality of 
template-based p hoto -zs measured using the Coadd pho¬ 
tometry ('Section l5.4|l . Compared against spec- 2 s from 
VVDS, template-based redshift estimators (see Section 
applied to the Coadd delivers photo-zs of compara¬ 
ble quality + z) ss 0.05 with an o utlier fraction of 

^10% ) to those reported for CFHTLenS (iHevmans et al.l 
l2012f) . This is despite the fact that the PSF-homogenized 
CFHTLenS aperture photometry reaches more than a 
magnitude deeper than the Coadd profile-fit photome¬ 
try. 

2 . 2 . UKIDSS Large Area Survey 

http://skyserver.sdss.org 


The Large Ar ea Survey (LAS) comp onent of UKIDSS 
is described by [Lawrence et al.l (|2007f) and utilizes the 
Wide Field Camera (WFCAM) on the 3.8 m United 
Kingdom Infrared Telescope (UKIRT). The LAS is tak¬ 
ing advantage of the large WFCAM field of view (0.21 
deg^) to image 4000 deg^ in the YJHK filter set. As 
of UKIDSS Data Release 10, the LAS was roughly two- 
thirds complete in terms of area covered. The depth goals 
in AB magnitudes are Y = 20.9, J = 20.4, H = 20.0, 
and K = 2 0.1. We discuss UKIDSS depths further in 
Section 17.11 

To construct the S82-MGC, we obtained public 
UKIDSS-LAS Data Release 8 (DR 8 ) catalogs from the 
WFCAM Science Archived (WSA). The LAS covers the 
majority of Stripe 82 and was substantially complete in 
this region in DR 8 , containing 5,072,574 near-IR selected 
sources, although some holes remained. 

When aperture magnitudes are required, we use 
UKIDSS magnitudes computed in the 4*^ aperture op¬ 
tion (AperMag4) with a radius of 2'.'8, unless otherwise 
noted. Aperture magnitudes in UKIDSS catalogs have 
been adjusted to “total magnitudes” under the assump¬ 
tion that every object is a point-source. To handle ex¬ 
tended sources, a complex query (see Appendix) is re¬ 
quired to obtain the “correction” (e.g., AperCor4) so that 
it can be removed. We also apply Calactic extinction cor¬ 
rections by subtracting the UKIDSS tabulated AX values, 
(where X can be Y JHK). 

The WFCAM focal plane consists of four separate de¬ 
tectors arranged in a checkerboard pattern. UKIDSS ob¬ 
tains large survey coverage by tiling this pattern across 
the sky. In any given position, multiple short exposures 
are coadded into a data unit referred to as a multiframe. 

2.2.1. UKIDSS imaging depth and seeing 

The UKIDSS imaging depth varies with location in 
different ways for different bands. We use the LAS cata¬ 
log information to estimate the depth of the imaging in 
which each source was detected and measured as follows. 
First we identify the set of sources in a given band be¬ 
longing to a specified multiframe. We then compare the 
AperMag magnitudes against their magnitude errors and 
interpolate to find the magnitude at which the average 
error is 0.1 mag. We consider this a rough estimate of 
the 10(7 depth of that particular multiframe and asso¬ 
ciate this depth limit with all sources belonging to that 
multiframe. Some multiframes do not have enough de¬ 
tections to provide a solid depth estimate; this situation 
is often an indication of a problem in the UKIDSS data. 
Sources from such multiframes are assigned a magnitude 
depth of 0.0. The depth information is used to define a 
flux-limited survey area in Section 0 

We obtain the seeing measurement (FWHM) in each 
UKIDSS band from the MultiframeDetector tables on 
the WSA. See the Appendix for the relevant query. 

2.3. SDSS Spectroscopic Redshifts 

Our primary source of spectro scopic redshifts is from 
BOSS, an SDSS-HI program (jEisenstein et al.l 1201111 
that obtained spectroscopic redshifts for 1.5 million 
galaxies over 10,000 deg^ using the BOSS spectro- 

http://surveys.roe.ac.uk/wsa 













































4 


Bundy et al. 


graphs (|Smee et al.ll2012D on the 2.5 m Sl oan Founda¬ 
tion T elescope at Apache Point Observatory (|Gunn et al.l 
200611. An overview of the BOSS survey is presented in 
Dawson et al.l (j2013[l . BOSS gal axies were selected from 
SDSS Data Release 8 imaging (jAihara et al.ll2011^ us¬ 
ing a series of color and magnitude cuts. A concise re¬ 
view of the selection criteria for the two BOSS samples 
is presented in Paper II. The LOWZ sample targeted lu¬ 
minous red galaxies (LRGs) at 0.15 < 2 < 0.43, while 
the GMASS sample aimed to collect massive galaxies 
in a broader color range at 0.43 < z < 0.7. Redshifts 
were measured usin g the processing pipeline described in 
iBolton et al.l (|2012ll . We also include redshifts of SDSS- 
I/II “Legacy” objects, primarily LRGs. 

All SDSS redshifts in the S82-MGC were obtained by 
cross-matching to the SpecObj-drlO cataloeFI using a O'.'S 
tolerance. Matching is done only to specObj entries that 
have SPEGPRIMARY=1 and ZWARNING.NOQSO=0 
(for BOSS programNcunes) and ZWARNING=0 (for non- 
BOSS programNcunes). Both stars and galaxies are 
matched, producing a total of 149,439 spectroscopic red¬ 
shifts included in the S82-MGC. 

3. PHOTOMETRY 

3.1. Matched ugrizY JHK SYNMAG Photometry 

Photometry matched to a consistent point-spread func¬ 
tion (PSF) across the SDSS-I-UKIDSS filterbands is a 
key component of the S82-MGC that enables SED fitting 
and near-IR based M, estimates. We describe our use of 
catalog-level matc hing SYNMAGs in thi s section. 

As described in iBundv et al.l (j2012t) , consistent pho¬ 
tometry across multiple filter bands with varying PSFs 
(and especially from multiple telescopes) is most often 
obtained by convolving the images available in each band 
to the poorest PSF in the ensemble. The same aperture 
is then applied to a source detection in every image and 
a consistent measure of the light from the same regions 
of extended objects at different wavelengths is obtained. 
While robust, this method reduces the image quality of 
all data sets and requires substantial effort in under¬ 
standing and processing the imaging pixel data. The 
task becomes more challenging when matching data sets 
over large areas as data volumes extend to the terabyte 
regime. More sophisticated model “forced” fitting photo¬ 
metric techniques may provide better estimates but are 
even more difficult to execute over large areas. 

In an era of rapidly expanding and overlapping imag¬ 
ing surveys, a faster alternative that operates at the cat¬ 
alog level (instead of the pixel data) is clearly appeal- 
ing. This is the m otivation for the SYNMAG technique 
(|Bundv et al.ll201§l that we apply here. 

Briefly, starting with the available SDSS and UKIDSS 
catalogs, we cross-match sources usi ng a position toler- 
ance of 0'.'8 and “handshaking” ('e.g.. lHewett et al1l2006f) 
to reduce misidentifications. Handshaking requires that 
the nearest UKIDSS source matched to an SDSS object 
have as its nearest match that same SDSS source. Thus, 
cross-referenced pairs of sources from the two catalogs 
are uniquely matched. From the original Goadd par¬ 
ent catalog with 15,342,585 sources, 3,175,036 (20%) are 

Available at http://www. sdss3.org/drlO/tutorials/lss_galaxy. 
No additional spectra were made available in the relevant 
fields in DR12. 


matched to a UKIDSS source. Of the original UKIDSS 
catalog, matched sources account for 63%. 

In the next step, SYNMAGs require that intrinsic, PSF- 
convolved surface brightness profiles have been fit in at 
least one band (we will refer to it as the “profile band”). 
In our case, we use the r-band profile fits measured by the 
SDSS pipeline, which also delivers the SDSS ModelMags. 
Given the seeing and photometric aperture in a target 
band, we use the profile fit to “predict” the PSF-matched 
aperture photometry that would have been observed in 
the profile band. This approach provides a consistent 
color between these two bands. We then repeat the 
process for the remaining profile-target band combina¬ 
tions, thereby building a full set of matched photome¬ 
try. By initially decomposing the profile fits into Gaus¬ 
sian mixture models (also see lHogg fc Landl2013l l. PSF- 
convolutions can be done analytically, and because all 
required information is derived from source catalogs, the 
technique is extremely fast. Further detail s and tests of 
the methodology are given in iBundv et al.l (|2012h . 

3.2. UKIDSS Total Magnitudes for Extended Sources 

Measurements of total magnitudes for extended 
sources are problematic in the UKIDSS public catalogs. 
The WSA “known issues” pag43 reports several long¬ 
stand ing problems. Kron radii and magnitudes (|Kronl 
Il980f) are not reliable in any survey data except for the 
Ultra Deep Survey. Other total flux estimators such 
as Pe trosian (jPetrosianI 119761 1 and Hall (jHall fc Mackavl 
119841 1 fluxes, as well as their associated radii, are not 
measured correctly for blended sources, as we demon¬ 
strate below. Unfortunately, while the overall fraction 
of UKIDSS sources that are affected by blends is small 
(about 1% among all AT-band detections in the S82- 
MGC), the fraction rises to 8-10% for galaxy samples 
with logM*/M 0 > 11.4 and is magnitude (i.e., redshift) 
dependent. For blended sources in the S82-MGC, the so¬ 
lution we adopt is to apply a total magnitude aperture 
correction to each UKIDSS band that is equal to the 
difference between the synthetic aperture z-band pho¬ 
tometry and the z-band CModelMag total flux estimator. 
In other words, we will assume that the total flux aper¬ 
ture correction in the z-band for these objects applies 
to the near-IR bands as well. Similar assumptions have 
been tested and a pplied in other surveys as well (e.g., 
iCaoak et al.ll2007^ . 

We begin by demonstrating the amplitude of the 
UKIDSS total magnitude errors. We have downloaded 
nearly 6000 UKIDSS image “cutouts” centered on K- 
band bright sources, both isolated and blended and 
drawn from the S82-MG C, each with a size of 1' x 1'. 
We have run SExtractor (iBertin &: Arnoutill996ll on the 
cutouts and compared its output magnitudes to various 
total flux estimators in the UKIDSS DR 8 catalog. The 
results are reported in Figure [TJ The zeropoint is deter¬ 
mined from the keywords stored in the zeroth and first 
extension hea ders in each cuto ut FITS image, following 
Equation 1 of lHill et al.l (|2011l l: 

ZP = ZPmag + 2.51ogf —ExtX ((AMi-|-AM 2)/2 —1) (1) 
where ZP^ag refers to the keyword, MAGZPT, t refers 
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Fig. 1.— Comparisons of various total magnitude estimates from the UKIDSS LAS catalog to MAG_AUTO remeasured with SExtractor. 
UKIDSS Petrosian magnitudes (Petro), Kron magnitudes (Kron), isophotal magnitudes (Iso), Hall magnitudes (Hall), and corrected “total” 
AperMagd magnitudes (Aper6) are displayed. Hatched histograms refer to offsets for isolated sources while solid, filled histographs show 
results for deblended sources. 


to EXP.TIME, Ext refers to EXTINCT, and AMi and 
AM2 refer to AMSTART and AMEND, the beginning 
and ending airmass values of the exposure. We verify 
our zeropoint determination by first comparing aperture 
photometry in the UKIDSS AperMag2 through AperMagS 
apertures (corresponding to aperture radii of 1.4, 2, 2.8, 
4, and 5.7 arcseconds). Focusing on isolated sources 
and de-correcting the LAS aperture magnitudes as in 
Section 13.11 we find that the median LAS photometry 
is '^0.03 magnitudes brighter in all apertures; we at- 
tri bute this to small ca talog-level corrections discussed 
in iHodgkin et al.l (|2009ll . 

Using hatched histograms for isolated sources and 
solid, filled histograms for blended sources. Figured] dis¬ 
plays blended total photometry errors in the UKIDSS 
Petrosian magnitudes (Petro), Kron magnitudes (Kron), 
Hall magnitudes (Hall), and corrected “total” AperMagS 
magnitudes (Aper6). In all cases except for isopho¬ 
tal and AperMagS magnitudes, LAS catalog photome¬ 
try of deblended sources is significantly brighter than 
MAC_AUTO, typically by 0.5-1 magnitude. Petrosian 
magnitudes provide the poorest estimates, with offsets 
as large as 3 magnitudes. 


For isolated sources, Hall magnitudes correlate best 
with MAG_AUTO but still suffer biases in the case of 
blended sources. The “total” magnitudes reported for 
AperMagS (i.e., with the UKIDSS aperture correction 
applied) show the opposite behavior in that blends are 
fainter than MAG_AUTO on average compared to iso¬ 
lated sources. Even when comparing de-corrected AperS 
magnitudes to the corresponding aperture photometry 
remeasured by SExtractor, there is a bias such that the 
reported AperMagS photometry is 0.05 mag fainter for 
blends. This bias decreases as the aperture decreases 
and nearly disappears for Aper4 and smaller apertures. 
We speculate that the reported aperture magnitudes in 
the UKIDSS catalog have also been adjusted to remove 
an estimated contribution from overlapping sources in 
blended objects, although we could not find UKIDSS 
documentation to confirm this conclusion. The compari¬ 
son to “total” AperMagb values suggests this adjustment 
introduces a roughly 0.1 magnitude bias towards fainter 
magnitudes in the large-aperture photometry of blended 
sources. 

The offset for isophotal magnitudes shows no appar¬ 
ent dependence on whether a source is blended. A 
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Fig. 2. — Illustration of the UKIDSS total magnitude estimator 
adopted in the S82-MGC. For sources classified as galaxies we plot 
it’-band results in this figure, but the other bands are similar. Here 
CModelMag refers to the UKIDSS aperture magnitude corrected by 
the difference between the z-band synthetic aperture and the z- 
band CModelMag. Blended sources are highlighted in purple and 
fall off sharply beyond the detection limit (~19.5). In the top 
panel, the Hall magnitudes are fainter than CModelMag for bright 
galaxies. The green curve is a fit to this trend. In the bottom 
panel, we have corrected the Hall mags of isolated sources using 
this fit and substituted the CModelMag value for blended sources. 
This plot of the Hall mag as the dependent variable confirms that 
there are no biases in the total flux estimator above the detection 
limit. 

simple brightening of the reported isophotal photom¬ 
etry by ~ 0.2 magnitudes brings them into agreement 
with MAG_AUTO, but the relatively large scatter in 
the isophotal offset hints at the obvious problem that 
isophotal magnitudes cannot account for surface bright¬ 
ness dimming and will introduce biases as a function of 
magnitude (or redshift). 

Without a robust total magnitude estimate from 
UKIDSS for blended sources, we turn to SDSS. Not only 
does the SDSS Coadd photometry deal with blends, it 
also fits 2D surface brightness profiles to every source. 

A recommendecEB total flux estimator is the CModelMag 

see http://www.sdss.org/drl2/algorithms/magnitudes/#which 


which takes the best linear combination of separate 
de Vaucouleurs and exponential fits in every SDSS 
band. Such profile fitting can account for extended 
surface brightness profiles that cross under the noise 
background, but, constrained by the profile shape as- 
sumptions adopted, they can also introduce biases (e.g., 
iBernardi et al.ll20l^ . 

With these caveats in mind, we build a new total 
magnitude estimate, HallTot, that is referenced to the 
UKIDSS Hall magnitudes for non-blended sources and 
to the SDSS z-band CModelMag magnitude for blends. 
For blends, we apply the offset between the z-band syn¬ 
thetic aperture magnitude and CModelMag^ to the aper¬ 
ture magnitude (AperMag4) measured for each UKIDSS 
band. This choice anchors our UKIDSS total magnitude 
estimates to CModelMag^. However, we would prefer to 
use the UKIDSS Hall mags for isolated sources, given 
their strong performance with respect to the remeasured 
MAG_AUTO. The solution is to adjust the Hall mags to 
also match CModelMag^ on average. The process is illus¬ 
trated in Figure [2] and delivers blend-resistant HallTot 
magnitudes for all four UKIDSS bands. 

4. STAR-GALAXY SEPARATION 

Star-galaxy separation in the S82-MGC is more com¬ 
plicated than for typical SDSS imaging data because 
of problems in the PSF characteri zation of the Coadd 
(|Annis et al.ll20l4 iHuff et al.ll2014fl . Instead of a classi¬ 
fication based on the SDSS TYPE parameter, we use SDSS 
shape information as a crude first estimate and then re- 
hne the classification using the (J — K) vs. (5 — i) colors 
to separate g alaxies (and q ua sars) from th e stellar locus. 

Following iBaldrv et al.l (I2010D and iStrauss et al.l 
( 1200 ^ , we define the SDSS star-galaxy separation pa¬ 
rameter as Asg = Tpsf - Tmodei, where rpsf and rmodei are 
the r-band PSF and ModelMag measurements, respec¬ 
tively. Sources with Agg > 8 can be confidently classified 
as (large) galaxies and those with Agg < 0.05 as point 
sources. For the remaining sources, we adopt a modified 
version of the (J — K] distan ce from the stellar locus 
defined by IBaldrv et al.l (I 2010 D : 


Aggjk — Ar)AB /locus(l/ 0 (^) 

where we correct all colors for Galactic extinction and 
/iocus (5 -i) = -0.523 for {g - i) < 0.7, /iocus (3 - *) = 
—0.1632 for {g — i) > 2.3. For 0.7 < {g — i) < 2.3 the 
locus is defined as: 

fiocus{g-i) = -0.89-P0.615*(g-i)-F0.13 +( 5 -i)^ (3) 

We adopt a color cut defined as Aggjk = 0.25 which is 
displayed in Figure |3l 

A check on the completeness and purity of the S82- 
MGC galaxy classification can be made by comparing 
to the star-galaxy separator (Leauthaud et ah, private 
communication) based on the overlapping GFHT CS82 i- 
band imaging survey (Kneib et ah, in preparation) which 
delivers excellent shape-based classifications thanks to a 
median seeing of FWHM=0'.'6. First, we note that dif¬ 
ferent internal classifications using CS82 are impure or 
incomplete at the several percent level, representing a 
quality floor in typical ground-based data sets. Next, we 
-rnaggympare against the “DW” GS82 classification which is 
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Fig. 3. — Top panel: Color-color criteria for star-galaxy sepa¬ 
ration adopted in the S82-MGC. Shaded contours and data points 
represent the distribution of sources with i < 21.5. The colors are 
Galactic extinction corrected. Bottom panel: Same as above but 
for SDSS Coadd TYPE=3 sources classified as galaxies. The TYPE 
based galaxy classification is contaminated by a significant fraction 
(~10%) of stars, likely as a result of PSF characterization prob¬ 
lems in the Coadd. SDSS Coadd point-sources (TYPE=6) do not 
contaminate the color-color based galaxy locus. 

based on manual divisions made in the size-magnitude 
plane. The completeness of the S82-MGC galaxy classifi¬ 
cation compared to DW is 96%, with a 1% level of con¬ 
tamination from stars. These numbers may be optimistic 
because both classifiers can fail. A visual inspection of 
SDSS images suggests that the s82-mgc galaxy sample 
is contaminated by stars (especially binaries) at the few 
percent level. 

While our classification requires a (J — K) color, it 
performs much bette r than the Coadd TYPE parameter, 
used, for example, bv iReis et al.l (|2012ll to define a galaxy 
photo-z sample. As shown in the bottom panel of Fig¬ 
ure [3l a Coadd selection with TYPE=3 leads to a galaxy 
sample with at least 10% contamination from the stellar 


locus. 


5. PHOTOMETRIC REDSHIFTS 

Despite the large number of spectroscopic redshifts in 
S82-MGC (nearly 73,000 galaxy spec-zs, or 11% of galax¬ 
ies with i < 22.5), reliable photometric redshifts are 
paramount for building a complete sample in Stripe 82, 
given that the majority of the spec-zs were obtained from 
SDSS campaigns employing complex selection criteria. 

We have explored a number of photometric redshift 
estimators and present plots summarizing their perfor¬ 
mance against available spectroscopic redshifts in Fig¬ 
ures [4] and [3 We break with tradition and characterize 
photo-z scatter without dividing by (1 -I- z) since there 
is no evidence for such trends in our data. We quote 
(7z as the 3(T-clipped standard deviation and characterize 
the outlier fractions with separate parameters. We define 
catastrophic outliers as those with |Az| >0.1, where Az 
is difference between the photometric and spectroscopic 
redshift. Statistical outliers refer to photo-zs with |Az| 
greater than ScTz- 

The best photo-z performance for massive galaxies 
in clusters is produced by the red-seque nce Matched- 
filter Probabilistic Percolation (redMaFFer. IRvkoB: et all 
[MD. For non-cluster galaxies, the best performance 
is produced by the r ed-sequence Matc hed filter Calaxy 
Catalog {redMaGiC, iRozo et ^1201511 and neural net- 
work results derived in iReis et al.l ( 20121) using ANNz 
(|Collister fc Laha'^l2004ll . We use a combination of these 
estimates to produce the Zbest redshift estim ator appro¬ 
priate for massive galaxies (see Section [0|l . 

We also coni puted templat e-based photo-zs using BPZ 
(|Benftej 1200011 and EAZY (|Brammer et al.l 1200^ (see 
Section 15.41) . Even after extensive experimentation with 
various template sets, priors, and zero point offsets, we 
could not achieve results that were competitive with 
ANNz (let alone redMaFFer) for bright galaxies. The 
scatter (ctz) in photo-z-spec-z residuals was typically 
twice as large with the template codes as compared to 
ANNz, and the fraction of catastrophic outliers higher 
by roughly 50%. However, for faint galaxies {i > 22.5), 
template codes perform somewhat better, likely because 
of the smaller training sets available in the face of poorer 
photometry. We use template photo-zs as a check on the 
neural-network trained redshifts and additionally include 
them in the S82-MGC with a description below for stud¬ 
ies that can benefit from any redshift constraints for the 
faintest objects in the catalog. 


5.1. redMaFFer and redMaGiC photo-zs 

The beslE3 photo-zs we have compared—with CTz and 
catastrophic outlier fractions at the 1% level—are those 
from the recent redMaFFer and redMaGiC catalogs. 
However, these are limited to red cluster galaxies only 
{redMaFFer) or red g alaxies in general ( redMaGiC). As 
described in detail in iRvkoff et al.1 (j201-^ . redMaFFer’’s 
primary goal is to use the red-sequence method to ro¬ 
bustly identify galaxy clusters and their richness using 
imaging data sets. By training on available spec-zs, 
redMaFFer iteratively defines a model for the redshift 

A complete and fair comparison against ANNz would con- 
sider just the (primarily red) galaxies for which redMaPPer and 
redMaGiC photo- 2 :s are considered trustworthy. 
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Fig. 4.— Comparisons of three photometric redshift estimators to available spectroscopic redshifts in Stripe 82. The comparison is 
limited to z < 22.5 and 0.01 < 2spec < 0.8. The left and middle panels ar e from the redMaPPer project (Rozo et al., in preparation), while 
the right panel compares neural-network photo-zs from [Reis et al.l (I2012lh The 3cr-clipped dispersion is listed in each panel along with the 
fraction of catastrophic outliers defined by \Az\ > 0.1. Contours are plotted at high data densities with 0.3 dex logarithmic spacing in the 
left and middle panel and 0.4 dex in the right panel. The 1-to-l relation is plotted in each panel as a thin light grey line. 



Fig. 5.— Photometric redshift-dependent photo -2 scatter and 
outlier fractions for the photo -2 estimators used in this paper as 
compared to available spec- 2 s in Stripe 82 with i < 22.5 and 0.01 < 
•Zspec < 0.8. The 3f7-clipped dispersion {(Tz) is indicated by the 
solid lines and symbols. The redshift error has not been divided 
by (1 -|- z). The dashed lines plot /outlier? the fraction of photo- 2 s 
that deviate by more than 3cr; this is a more stringent definition 
than /catastrophic Sensitive to systematic patterns of scatter that 
can bias certain redshift bins. 

dependent red-sequence that can be used to deliver pre¬ 
cise photo-zs for red galaxies. 

We compare redMaPPer and redMaGiC photo-zs to 
available spec-zs in the first two panels of Figure 4. 
redMaPPer photo-zs are only available for cluster galax¬ 
ies, which we define as galaxies in the redMaPPer clus¬ 
ter members catalog that have a minimum membership 
probability of 90%. For these galaxies, we assign a pho¬ 
tometric redshift equal to the photometric redshift of 
the redMaPPer galaxy cluster. Roughly 4% of S82-MGC 
galaxies with i < 22.5 fall in this category. The second 
estimate we refer to as Zred2 (or ZREd2, middle panel of 
Figure 4). Zi.ed 2 is the starting point for the redMaGiC 
photometric redshifts. It is the result of applying the 
empirically derived and redshift-dependent red-sequence 
template from the cluster finding algorithm to all red 
galaxies in the sample, regardless of cluster membership. 
Approximately 20% of S82-MGC galaxies with i < 22.5 


have realizable Zred 2 estimates. Here, we adopt Zred 2 
estimates if the quality of the redMaPPer template fit 
satisfies < 10 and the galaxy is brighter than Lf.{z) 
luminosity used by the redMaPPer and redMaGiC algo¬ 
rithms. Our criteria for selecting redMaPPer and red¬ 
MaGiC redshifts were chosen to balance the number of 
usable photo-zs against their quality. More liberal cuts 
only marginally increase the number of available redMaP¬ 
Per and redMaGiC photo-zs, but significantly decrease 
their quality. 

5.2. ANNz photo-zs from Reis et al. 

IReis et all 1)201211 co mpute photo-zs using th e ANNz 
neural network code (jCollister fc Laha\I 120041 1 applied 
to r < 24.5 sou rces from the SDSS Coadd catalog 
()Annis et al.ll20l3l with Coadd TYPE=3. As described 
in Section S) using this Coadd TYPE criteri£(3 to select 
galaxies results in a sample with ^10% contamination 
from stars. A more accurate star/galaxy separation is 
presented in Section 0] and used to define a complete 
sampl e, ukwide, i n Sec tion [T] 

The IReis et all (1201211 photo-zs (we refer to them as 
•ZReis) were trained and validated with ^83,000 spec-zs in 
Stripe 82. They perform well compared to BOSS galax¬ 
ies despite having been trained by a more limited spec-z 
sample that included only 6,682 BOSS redshifts (Fig¬ 
ure [H right-most panel). Comparing to all spec-zs now 
available, we find = 0.028 and a catastrophic outlier 
fr action of 5.7%, in line with the performance reported 
in IReis et all (|2012ll . 

While the training set included additional spec-zs from 
deeper surveys, including DEEP2 and VVDS, it is pos¬ 
sible that ZReis estimates perform poorly for bright, pri¬ 
marily “blue” galaxies that are not well represented in 
the training but could be an important contribution to 
the total massive galaxy population. We test this pos¬ 
sibility in Figure |6] by comparing ZReis photo-zs to esti- 
mat es fr om template codes, zeazy and zbpz (see Sec¬ 
tion 15.41) . The comparison is made for sources classi¬ 
fied as galaxies in the S82-MGC with XogM^jM^ > 11.3 

Compac t galaxies misclas sified as point sources would also be 
missed in the IReis et a nunii) photo -2 sample. 
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Fig. 6. — Independent check of neural network ANNz photo-zs, ^Reis? against EAZY photo-zs (left) and BPZ photo- 2 :s (right) for a galaxy 
sample with no other redshift information available (roughly 30% of the UKWiDE sample across all redshifts). Catastrophic problems in 
2^Reis resulting from inadequate training sets are not evident. Hints of potential photo- 2 : biases at the 0.1 level follow different patterns in 
the two panels and are difficult to confirm because, by definition, no spec- 2 :s are available for this selection. 


and belonging to the ukwide subsample described be¬ 
low (Section [7]). We only plot galaxies for which a spec- 
z or redMaPPer photo-z is not available, thus focusing 
where our best redshift estimate, Zbest, comes from ZReis- 
About 30% of UKWIDE galaxies fall in this category at 
these masses. 

Figure | 6 ] demonstrates that there are no catastrophic 
problems in ZReis resulting from inadequate training sets 
for massive galaxies with z < 0.7. While there are hints 
of potential photo-z biases at the 0.1 level, the patterns of 
these biases are different when comparing to zeazy (left 
panel) or zbpz (right panel). Civen that all three esti¬ 
mators may have biases, it is difficult to say more about 
biases in ZReis photo-zs at this level without additional 
spec-zs. 

5.3. Combined photo-zs for bright galaxies: Zbest 

For bright galaxies {i < 22.5), we combine spec-z and 
photo-z results together into a single “best” estimate 
that we refer to as Zbest- In order of priority, Zbest is 
set to: 

1. Zspec: Spectroscopic redshifts passing quality cuts 
from SDSS (including “legacy” and BOSS sam¬ 
ples), VVDS, or DEEP2. 

2. zrm: redMaPPer photo-z for cluster members. 

3. Zredz: redMaCiC photo-z for red field galaxies. 

4. ZReis: ANNz photo-zs from iReis et al.l (|2012ll . 

The fractional contribution of these estimators as a 
function of both M* and redshift is shown for the ukwide 
sample in Figure [7l Paper II studies the SDSS-I-BOSS 
spectroscopic completeness evident in this figure. Stud¬ 
ies using SDSS-t-BOSS samples to z ~ 0.7 are more 
than 80% spectroscopically complete at logM*/M 0 > 
11.6 (right-most panel) and only moderately worse at 
logM*/M 0 >11.4 (middle panel). 


5.4. Template photo-zs for faint galaxies 

Templa te photo- zs from EAZY (|Brammer et al.l[2008[l 
and BPZ (|Benitezll2000ll have also been computed to pro¬ 
vide a check on ZReis (see above) and additional redshift 
informatio n for faint sources in the S82-MGC {i > 22.5). 
Unlike the iReis et ^ (j2012f l ANNz catalog, we provide 
template photo-zs regardless of star/galaxy classffica- 
tion, in part to enable studies of possibly misclassified 
compact galaxies, but also to allow future improvements 
in separating stars and galaxies. As we show, the ZReis, 
-Zeazy, and zbpz estimates fail in different ways at the 
faint end of the S82-MGC sample. One advantage of both 
EAZY and BPZ is that some control of the failure rate is 
possible with cuts on the output quality ODDS parameter, 
if one is willing to limit the usable sample size. These 
features make the template photo-zs useful for source 
galaxies in weak gravitational lensing studies, for exam¬ 
ple. 

It is well known that template errors in the rest-frame 
near-IR can lead to degraded ph oto-z performance whe n 
near-IR data is in cluded (e.g., iBrammer et al.l l2008f) . 

I Bundy et al.l (l2012f) demonstrated acceptable photo-z 
performance with near-IR SYNMAGs included as a way of 
validating the SYNMAG methodology. We confirm that 
the addition of UKIDSS photometry here does not im¬ 
prove the template photo-zs, but also does not degrade 
them when the photometry is appropriately weighted to 
take template errors into account. In what follows, both 
-Zeazy and zbpz estimates are therefore based solely on 
the optical ugriz SDSS Coadd (ModelMag) photometry. 

For the zeazy estimates, we built a custom grid of 
redshift priors as a function of apparent magnitude via 
comparisons to the SDSS-I-BOSS spec-zs at bright mag¬ 
nitudes and the complete VVDS spec-zs at faint magni¬ 
tudes. We converted SDSS Coadd inverse h yperbolic sine 
(asin h) magnitudes (or “luptitudes”, see iLunton et al.l 
Il999tl to fluxes and used these as input without any ad¬ 
justments to the filter band zeropoints. We found the 
best performance at faint magnitudes using the standard 
set of EAZY templates. For bright galaxies (i ^ 20)— 
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Fig. 7.— Fractional contribution of different redshift results as a function of M* and redshift to the “best” estimator, 2 :best 5 appropriate 
for bright galaxies (i S. 22.5) in the S82 -MGC. The red line (zspec) indicates spec -2 from SDSS. The labels 2rm and z^ed refer to redMaPPer. 
ANNz photo- 2 S from IReis et al.l II2Q12I) are labeled 2phot ^rid shown in green. See text for details. This figure is reproduced and described 
further in Paper II. 


the relevant sample for checking the ANNz ZReis redshifts 
in Figure E] —we found be tter results with variants of the 
iBlanton fc Roweii (|2007f) templates (BR07). Aiming to 
improve outliers in the troublesome Zgpec ~ 0.35 range, 
we combined results from the EAZY-supplied BR07 tem¬ 
plate set with a modified version that removed the most 
extreme star-forming template. We further attempted 
to calibrate any remaining zeazy biases and refer to this 
version of zeazy estimates as “norm_brt”. Compared 
to SDSS-I-BOSS spec-zs with no ODDS cut, the norm_brt 
sample exhibited a Uz of 0.03 with an outlier fraction of 
18% and a catastrophic outlier fraction of ^10%. 

To estimate zbpz, no adjustment to the assumed pri¬ 
ors was made, but we did apply the zeropoint offsets 
that BPZ is able to compute for sources with spec-zs. 
We experimented with several different magnitude cuts 
when determining these offsets and also tested several 
SED template sets. These experiments led to only small 
improvements in photo-z performance. 

Given the superior performance of the redMaPPer and 
ANNz redshifts described above for bright galaxies, we 
turn now to the quality of zeazy and zbpz estimates 
for the faintest galaxies detected in the S82-MGC. Here 
z^EAZY is derived using the standard EAZY template set. 
Figure [5] presents comparisons to VVDS and DEEP2 
spec-zs. VVDS is a magnitude limited survey with 
i < 22.5, while DEEP2 targets were color-selected to 
have z > 0.7 and typically have magnitudes of i > 22.5 
which therefore reach or exceed the S82-MGC detection 
limits. Catastrophic outliers in the zeazy estimates tend 
to congregate near zeazy'^ 1-1, likely reflecting the cus¬ 
tom priors we have adopted. The zbpz outliers are more 
distributed and tend to fall below the true redshift. Par¬ 
ticular science cases may prefer one or the other failure 
mode. In addition, cuts on the ODDS photo-z quality pa¬ 
rameter allows cleaner photo-z samples to be obtained 
(lower panels). For comparison, the ZReis estimates for 
the DEEP2 sample are shown in Figure IHl 

6. SURVEY GEOMETRY 

A precise measurement of the total valid solid angle 
over which the S82-MGC sample was drawn is obviously 
critical to any derivation of source densities. Unfortu¬ 


nately, this measurement is complicated by the relatively 
large number of intersecting data sets: the Stripe 82 
Coadd imaging, the UKIDSS photometry in four sep¬ 
arate bands, and the BOSS spectroscopy. Not only must 
we define the geometric footprint of each data set but 
also the regions that should be rejected for various rea¬ 
sons, including image artificats and bright stars. 

To do so, w e have relied heavily on the MANGLE soft¬ 
ware package (| Swanson et al.l[200M^ supplemented with 
both custom built routines and IDL procedures bundled 
with the SDSS idlutil^^ package. With mangle it is 
possible to create, manipulate, and determine the inter¬ 
sections and overlap of spherical polygons on the sky. 
We refer to a set of such polygons as a “mask,” which 
can represent a survey’s on-sky footprint or a series of 
“holes” where data should not be considered. Due to 
the presence of a bug in combining polygon weights in 
previous MANGLE versions, we used a development ver¬ 
sion downloaded in June 2014 from the MANGLE GitHub 

repositorj|3- 

All MANGLE polygon files in S82-MGC are “pixelized” 
to a fixed resolution of r = 11 and “snapped.” We 
adopt angle tolerances following the DES collaboration 
(Aurelien Benoit-Levy, private communication) with cap 
tolerance variables, “a,” “b,” and “t,” set to 070027 and 
the intersection angle tolerance, “m,” set to 10“®". From 
a visual inspection of output polygon masks, we estimate 
that these mangle settings yield total area estimates 
that are accurate at the 2% level. A small ~0.8 deg^ 
portion of the final polygon file, with source detections 
overplotted, is presented in Figure ITUl 

6.1. BOSS Survey Masks 

Because the BOSS team has used mangle tools ex¬ 
tensively, several useful masks defining the BOSS sur¬ 
vey already exist and are publicly available. We used 
masks described on the BOSS DR9 Large Scale Struc¬ 
ture pag43 on the SDSS-III website. There is little to 

http://space.mit.edu/molly/mangle 
http://www.sdss.org/drl2/software/idlutils 
https://github.com/mollyswanson/mangle 
https://www.sdss3.org/dr9/tutorials/lss_galaxy.php 
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Fig. 8. — Photometric redshift quality for the faintest galaxies in the S82-MGC. Results using the template-based EAZY (left panels) 
and BPZ (right panels) photo -2 codes are compared to VVDS and DEEP2 spec- 2 s as indicated. VVDS is a magnitude limited survey 
with i < 22.5, while DEEP2 targets were color-selected to have 2 > 0.7 and typically have magnitudes of i > 22.5 which therefore reach 
or exceed the S82-MGC detection limits. Catastrophic outliers display significantly different behavior between the two codes. Given sparse 
training sets at faint magnitudes, the ODDS parameter output from EAZY and BPZ provides an advantage over neural-network photo- 2 s 
for selecting better performing photo -2 samples. For example, the lower panels show improved results for 0DDS>0.7. The percentage of the 
original sample that satisfies the ODDS cut is labeled in each panel as are metrics characterizing the resulting photo -2 quality. 


no difference in the BOSS DRIO masks because no new 
observations were obtained in Stripe 82 after DR9. Be¬ 
cause the BOSS footprint overlaps almost entirely with 
S82-MGC, applying the BOSS masks when selecting the 
final S82-MGC sample enables studies of the complete¬ 
ness of BOSS samples (Paper II) and reduces the total 
area by less than 10%, almost entirely as a result of the 
BOSS collision priority rejection mask (see Section lOl) . 

The BOSS survey footprint (acceptance 
mask) is defined in the mangle polygon hie, 
boss_geometry_2011_06_10 .ply. We also make 
use of several rejection masks. Beyond the bright star 
rejection mask, bright_star_mask_pix.ply, the cen- 
terpost mask, centerpostjnask.ply, removes a small 


region associated with the center of each BOSS plate 
where spectroscopic targets could not be observed. The 
collision priority mask, collision_priorityjnask.ply, 
accounts for small zones around BOSS galaxy targets 
that may have been de-prioritized in favor of other 
source classes (such as quasars). Before application to 
S82-MGC, all BOSS masks are trimmed to the Stripe 82 
region. 

6.2. UKIDSS footprint masks 

We derive the areal geometry of the four hlter bands 
in the UKIDSS LAS photometry in Stripe 82 using cus¬ 
tom software and make the results publicly available on 
the S82-MGC website. We hrst dehne a template for the 
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Fig. 9.— Performance of the 2Reis neural-network photo- 2 :s when 
compared to the VVDS and DEEP2 samples used in Figure [8] 
The 2Reis estimates were trained in part on these spec- 2 s but show 
larger scatter (c^ than the template fitting photo- 2 : estimates com¬ 
pared in Figure[§] It is also more difficult to prune high-performing 
subsets of neural-network estimates. 


GRIDZONE_089, area = 0.8, frac = 0.82 
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Fig. 10.— A small portion of S82-MGC illustrating the union of 
acceptable survey polygons, each shaded with different colors for 
visual clarity. White regions correspond to rejection masks for 
stars, missing data, imaging artifacts, and BOSS targeting limita¬ 
tions. The positions of galaxies from the UKWiDE subsample (see 
Section [tJi are overplotted. 

4-detector WFCAM imaging footprint following the field 
layout information on the UKIRT websit<S, and then 
identify all catalog sources associated with a specific WF¬ 
CAM multiframe. Multiframes can overlap, miss data 
from certain detectors, or be truncated, for example from 
spilling over the edge of the SDSS Coadd region that de- 

http://www.jach.hawaii.edu/UKIRT/instruments/wfcam/user. 


fined the positional boundaries in the UKDISS catalog 
query. 

For these reasons, we use the catalog source posi¬ 
tions to reverse engineer each multiframe's sky posi¬ 
tion and determine the geometry of its acceptable imag¬ 
ing. We construct heavily smoothed histograms of the 
celestial coordinates of associated sources in order to 
identify the gaps that define the WFCAM footprint. 
We then adjust the footprint template to match, de¬ 
riving spherical rectangles that correspond to each de¬ 
tector. Where multiframes are truncated we use this 
gap information (if a gap exists) and the most ex¬ 
treme source positions to define the maximum extent 
of each detector. The resulting “framemask” is as¬ 
signed a “weight” equal to the derived AB magnitude 
imaging depth discussed in Section 12.2.11 and is con¬ 
verted into the MANGLE polygon format. The pixelized 
and snapped UKIDSS footprint masks have file names 
of the form, mangleJramiemask_{y, j ,h,k}_v???.ply, 
where v??? represents a version number. 

6.3. Rejection masks 

Even in regions with valid UKIDSS or SDSS Coadd 
photometry, it is important to mask regions around 
bright stars and associated diffraction spikes, satellite 
trails, and other image artifacts. Ideally, one would de¬ 
fine masks by inspecting actual imaging data, but for 
wide-field surveys with many bands, this is a challenging 
data management problem. Instead, we inspect the spa¬ 
tial distribution of cataloged sources across the survey re¬ 
gion, distinguishing valid photometric sources from those 
flagged because of a problem identified by the photomet¬ 
ric pipeline. Problematic regions can be identified eas¬ 
ily by eye because they exhibit an overdensity of “bad” 
sources or trace an unphysical geometric pattern (such 
as a diffraction spike). 

We developed a software package that takes a catalog 
of source positions as input and allows the user to scan 
through the survey footprint. Bad areas can be traced in¬ 
teractively. The user can draw either a circle by defining 
the center with one click and the radius with a second, or 
a generic polygon by clicking a set of points to define the 
polygon vertices. These shapes are saved, over-plotted, 
and can be deleted and redrawn. When finished, the soft¬ 
ware converts the stored shapes into a MANGLE polygon 
file that can be used as a rejection mask to define the 
final survey footprint. 

This technique was used to define the rejection masks 
for all four UKIDSS bands as well as the SDSS Coadd 
catalog. For each UKIDSS band, the amount of rejected 
area is ~3 deg^. Masking of the SDSS Coadd accounts 
for 1.2 deg^. Most bright stars were already masked by 
the BOSS bright star rejection mask (4.3 deg^ over s82- 
MGC), but several were missed or were more of a problem 
in the UKIDSS bands. In the S82-MGC region, the BOSS 
collision priority mask sums to 16.8 deg^ and the center- 
post rejection mask amounts to 0.1 deg^. The rejection 
masks are available on the S82-MGC website. 

7. THE UKWIDE SAMPLE 

With the parent SDSS Coadd and UKIDSS data sets 
defined above, we select a galaxy sample optimized to 
have a well characterized M* completeness limit as wide 
in the S82-MGC. We refer to this 
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sample as “ukwide.” Forming the basis for the num¬ 
ber density evolution and BOSS completeness studies 
described in Paper II, ukwide contains 517,714 galax¬ 
ies with matched Coadd-|-UKIDSS photometry. It cov¬ 
ers 139.4 deg^, and is complete above logM*/M0 Ri II.2 
at z = 0.7. Data tables corresponding to the ukwide 
sample are made publicly available as described in the 
Appendix. 

The UKWIDE sample applies several initial cuts. Galax¬ 
ies are separated from stars as described in Section 01 
We then remove sources located in any of the rejection 
masks, accounting for bad regions in all four bands of the 
UKIDSS photometry and the SDSS Coadd imaging, as 
well as bright stars and zones where BOSS spectroscopy 
was not possible or incomplete (see Section [6] for more 
details). We also require the ukwide galaxies to fall in 
the BOSS DR9 acceptance mask and pass cuts on the 
UKIDSS ppErrBits error flag. 

We must also define the acceptable UKIDSS imaging 
depths for the ukwide sample. In AB magnitudes, we 
choose 5(7 magnitude limits of [20.32 19.99,19.56,19.41] 
for the four bands, YJHK (see l2.2.f]) . ukwide contains 
only those sources drawn from UKIDSS imaging with 
depths greater than these values. The motivation for 
these limits is discussed further below. 

7.1. M* Completeness Considerations 

In constructing the ukwide sample, we have adopted 
a selection that maximizes the available area with rela¬ 
tively uniform completeness in stellar mass. A conserva¬ 
tive estimate for this redshift-dependent M, complete¬ 
ness limit in any filter can be obtained by considering 
the theoretical range of galaxy stellar mass-to-light ra¬ 
tios (M*/L) present at each redshift. If we choose a 
reasonable representative SED for the maximum M^/L 
possible (e.g., a stellar population that formed quickly 
in a burst of star formation at z ^ 5) and compare it 
against the limiting magnitude in the filter band, we can 
estimate the minimum M* of all galaxies at that redshift 
with brighter observed magnitudes. 

For a reference redshift, we choose z = 0.7, which 
roughly corresponds to the redshift limit of the BOSS 
spectroscopy and the point where the available photo¬ 
metric redshifts become increasingly imprecise. Figure 
[TT] displays the M* limits described above as a function 
of the magnitude in several bands at this redshift. We 
focus on the rizYJHK as they provide the most crit¬ 
ical constraints on M^/L and M^, at these redshifts in 
the S82-MGC data set. The ideal magnitude limits in 
each band at some desired M* limit are simply given by 
the intersection of the plotted lines and that M^, limit. 
The observed r-band limit is significantly fainter than 
the other bands, a result of the 4000A break at this red¬ 
shift and the red colors defining the extreme-M*/L SED 
template. 

The SDSS Coadd riz 90% completeness depths are 
indicated by the solid symbols on the corresponding 
riz lines. The r-band depth limits the completeness 
of any S82-MGC sample to log ps ll.l. The 

redder bands can afford brighter magnitude limits (in 
AB units) and still have deeper M* limits at z = 0.7. 
By roughly matching the M* limit defined by the r- 
band, we can define ideal magnitude limits (5 (t) for the 
UKIDSS bands. As shown again by the solid symbols for 



Fig. 11.— The estimated M* completeness limit as a function 
of the limiting depth in relevant S82-MGC filter bands at a specific 
redshift. Solid symbols indicate the locations of detection limits in 
each plotted band. The SDSS Coadd r-band depth is comparable 
to the UKIDSS depths in terms of M* completeness. 

the YJHK bands, our final adopted limits are indeed 
close to logM»/MQ « 11.1 at z = 0.7. In practice, the 
precise numbers ([20.32,19.99,19.56,19.41] for YJHK) 
were matched to a limit at somewhat lower redshift be¬ 
fore testing revealed the hdelity of S82-MGC to redshifts 
approaching 0.7. 

Because the UKIDSS imaging depth is not uniform 
across Stripe 82, a hnal consideration for the adopted 
magnitude limits is the corresponding amount of area 
available. Working in terms of the M* limit defined at the 
5(7 detection threshold, Figure[l2]displays the cumulative 
distribution of area imaged to depths greater than that 
limit for each UKIDSS band, again referenced to 2 = 0.7. 
A limit of log M^/Mq r: 11.1 at this redshift makes use 
of all the UKIDSS area available. 

We quote the final M* completeness limit for ukwide 
at logM*/M0 R 11.2. This is a conservative choice that 
acknowledges potential systematics in the way we esti¬ 
mate the UKIDSS and SDSS Coadd depths, as well as 
catastrophic errors that may occur in matched photome¬ 
try and SED fitting when using sources measured at the 
5(7 limit. Effectively, a limit of log /Mq r 11.2 corre¬ 
sponds to a brightening of our nominal magnitude limits 
by 0.25 magnitudes. 

8. DERIVED GALAXY PROPERTIES 
8.1. The birth parameter, feiooo 

The presence of recent star formation among galaxies 
in imaging data sets is often inferred from a rest-frame 
color that straddles the 4000 A break. A major bene¬ 
fit of near-IR data is the ability to use color-color di¬ 
agrams to distinguish red galaxies with underlying but 
dusty star formation from truly quiescent systems (e.g., 
IWilliams et al.l 12009(1 . Here we exploit the full optical 
through near-IR photometry in the S82-MGC to estimate 
a proxy for recent star formation that accounts for dust. 
We compress the 2D information in a color-color dia¬ 
gram into a single measurement by htting the full SEDs 
available in the S82-MGC with the KCORRECT package 
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Fig. 12.— We plot the area covered in each UKIDSS filter band 
as a function of the limiting depth expressed in terms of the M* 
completeness at a specific redshift. Assuming that all other bands 
are sufficiently deep, the completeness limit corresponds to the 
faintest theoretical galaxy with a maximum reasonable MjL that 
is detected at 5c7 in the relevant band. 
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Fig. 13.— Correlation between the 61000 birth parameter, a proxy 
for recent star formation, and position in a restframe optical-near- 
IR color-color diagram for a UKWiDE sample with 0.3 < ^best < 
0.6 and log M*/Mq > 11.2. The value of logbiooo is indicated 
by the associated color, as labeled. The most quiescent galaxies 
with the lowest 61000 parameters (colored red) are clustered at red 
{q — r) colors but rel atively blue (r — H) colors, as expected (e.g., 
IWilliams et al.l 120091) . Galaxies with more recent star formation 
(higher 61000 5 bluer colors) span a more extended range in this 
diagram. Dusty star-formers continue upwards into the optical red 
sequence, but can be distinguished by their redder (r — H) colors. 

(|Blanton fc Roweisll2007tl . KCORRECT provides an es¬ 
timate of the birth parameter 61000, defined as the ratio 
of the average star formation rate within the previous 1 
Gyr to the star formation rate averaged over the galaxy’s 
history. As with M* estimates, 61000 depends on the as¬ 
sumed models and priors used to the fit the observed 
SEDs. Figure [13] illustrates the correlation between 61000 
and location in the {g — r) versus (r — H) diagram. A 
tight grouping of galaxies with very low 61000 values is 


evident as is a more extended star-forming sequence that 
includes dusty systems that would otherwise be located 
on the (optical) red sequence. Visual inspection of the 
CS82 imaging confirms that these galaxies are predomi¬ 
nantly disk-like or disturbed. The use of a near-IR based 
61000 parameter allows them to be easily distinguished 
from truly quiescent systems. 

In the S82-MGC, galaxies at z ^ 0.55 and M* ^ 11.5 
have 61000 values as high as 0.7, suggesting an occasional 
high rate of recent star formation, but the vast major¬ 
ity are peaked near 61000 = 0 as expected for an old, 
passively evolving stellar population. 

8.2. Stellar Mass Estimates 

A key motivation of S82-MGC is providing near-IR pho¬ 
tometry for robust M* estimates since the range of pos¬ 
sible stellar M/L ratios for stellar populations decreases 
at near-IR wavelengths. Among the most massive galax¬ 
ies with z < 1, the expectation is that the majority have 
little star formation and are passively evolving. Thus, 
M* estimates b ased on optical ph otometry alone may be 
sufficient le.g.. iPforr et al.l[2012ll . Still, observed z-band 
corresponds to restframe g-band at z ^ 0.7, and even the 
observed z-band, which is typically shallower than r or 
i in single-epoch SDSS imaging, falls blueward of rest¬ 
frame r-band. Because young stars significantly bias M* 
estimates limited to blue restframe wavelengths, the fact 
that nearly 40% of the CMA SS sample is intrinsically as- 
sociated with the blue cloud (iMontero-Dorta et al.ll201^ 
and the evidence for an increase with redshift in the 
fraction of CMASS g alaxies with spect ra indicating re¬ 
cent star formation (jChen et al.l l2012| l motivate rest¬ 
frame near-IR photometry as an important ingredient 
for checking evolutionary resultS based on estimates of 
M*. 

8.2.1. S82-MGC fiducial M* estimates 

The fiducial M* estimates derived for the S82-MGC 
use th e Bayesian code initially pres e nted in iBundv et al.l 
(j2000[) and used in IBundv et al.l (j20ir)l l. Paper III 
presents additiona l mass estimates testin g various pri¬ 
ors with iSEDfit (iMoustakas et ahl 1201311 based on the 
S82-MGC photometry for the UKWIDE sample. The ob¬ 
served SED of ea ch galaxy is compared to a grid of 13,440 
models from the iBruzuai fc Charlotl (j2003[l (BC03) pop¬ 
ulation synthesis code that define a set of priors spanning 
a (randomized) range of metallicities, star formation his¬ 
tories (parameterized as exponentials), ages, and dust 
content. Ages are restricted to less than the cosmic age 
at ea ch redshift and no bursts are included. A Chabrier 
IMF (IChabrieill200^ . IIm= 0.3, nA=0.7, and a Hubble 
constant of 70 km s“^ Mpc“^ were adopted. 

At each grid point, the M^jL ratio in the reddest 
observed band, the value inferred from multiplying 
M^/L by the luminosity in this band, and the likelihood 
that the model matches the observed SED are stored. 
This likelihood is marginalized over the grid, giving an 

Paper III presents tests that, while revealing biases in optical 
M* estimates for some galaxies, show that the choice of stellar 
synthesis modeling and priors has a larger effect on evolutionary 
signals for log M*/Mq >11.3 and z < 0.7 than the use of optical- 
only M* estimates. 
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estimate of the stellar mass probability distributiorF^. 
We take the median as the final estimate of M*, which 
counts the “current” mass in stars and stellar remnants. 
The 68% width of the distribution provides an uncer¬ 
tainty value which is typically ^0.1 dex. 

8.2.2. Extant BOSS M* estimates 

Systematic uncertainties are commonly acknowledged 
in studies employing M* estimates but not always care¬ 
fully studied. The largest uncertainty is the IMF assump¬ 
tion which, to first order, can cause fa ctor of ^2 offsets 
i n values. Setting the IMF aside, IMoustakas et al.l 
(|20I3ll provide a careful examination at how different 
stellar population models and assumptions regarding pri¬ 
ors can influence M^, estimates based on th e same set of 
photometry. While IMoustakas et al.l (|2013ll demonstrate 
that the conclusions in their paper are largely insensitive 
to these systematic effects, we adopt a similar method¬ 
ology in Paper III that reveals the much greater impor¬ 
tance of adopted priors and assumptions for constraining 
massive galaxy evolution with high precision. 

This section presents a preview of Paper III by compar¬ 
ing the fiducial near-IR M^, estimates described above to 
publicly availablcp^ estimates for BOSS galaxies from the 
SDSS-III collaboration. For a more systematic study of 
why M* offsets can occur between estimators (including 
differences in the basis photometry set, stellar synthesis 
models, and adopted priors), please see Paper III. 

In the current work, we only make comparisons to pre¬ 
viously (spectroscopically) measured BOSS and SDSS 
Legacy galaxies that use either the single-epoch SDSS- 
only photometry or SDSS spectroscopy. We compare 
this d ataset to the “Wisconsin PCA” masses (jChen et al.l 
[Mil . tw o versions of th e “Por tsmouth” masses de¬ 
scribed in iMaraston et al.l (|2013l . referred to hereafter 
as M1 3), and the “Gra nada” masses released in SDSS 
DRIO (|Ahn et al.ll2014ll . The Granada and Portsmouth 
masses are based on stellar population synthesis SED 
fitting to the single-epoch SDSS ugriz ModelMag colors 
(after correcting for galactic extinction). The Wisconsin 
PGA masses are derived from a PGA analysis of opti¬ 
cal wavelength regions of the BOSS spectra. The re¬ 
sulting Mtf/L ratios in all three cases are scaled to to¬ 
tal magnitudes as given by the CModelMag_i measure¬ 
ments. All BOSS estimates assume a flat cosmology with 
Dm = 0.274 and Hq = 70 km s“^ Mpc“^. The impact on 
M* from using these cosmological parameters compared 
to those adopted for the S82-MGC M* estimates is less 
than 0.01 dex. 

The “passive” only version of the Portsmouth masses 
are derived by finding the best fit between the ugriz 
photometry and a single SED template built from the 
combination of two passively evolving, coeval bursts of 
star formation. The first burst is given solar metallicity 
and accounts for 97% of the total mass, while the second, 
accounting for the remainder, has 0.05 solar metallicity. 
Age is the only fitting parameter for this “passive” tem¬ 
plate but ages less than 3 Gyr are not allowed. A second 
set of M13 estimates, referred to as “SF”, are based upon 
fitting stellar population models with a range of star for- 

We assume the prior grid adequately samples the parameter 
space of the posterior. 

https://www.sdss3.org/drlO/spectro/galaxy.php 


mation histories, ages, and metallicities. Both exponen¬ 
tially declining and truncated star formation histories are 
included, and an upward correction of -1-0.25 dex is added 
based on an analysis of simulated data. Neither of the 
MI3 estimates includes possible dust extinction inside 
the target galaxy. We use DRIO versiontH with Kroupa 
IMFs for the MIS estimates. As with the S82-MGC es¬ 
timates described in Section 18. 2. II both M13 estimates 
measure the current mass in stars and ste llar remnants at 
the ag e of the stellar population model. IMaraston et al.l 
(|20I3ll combine the passive and SF estimates into a sin¬ 
gle, “best” estimate by choosing the passive estimate for 
galaxies with {g — i) > 2.35 and the SF estimate for 
{g — i) < 2.35. In comparisons described below, we at¬ 
tempt to reproduce the M13 “best” estimates by adopt¬ 
ing the passive or SF M* value from the Portsmouth 
catalogs using the Coadd {g — i) color. There will be 
slight differences from the actual M* values used in Ml3 
because the Coadd photometry is de eper. _ 

For the Wisconsin PCA masses, iChen et al.l (|2012[1 
compared their PCA eigenspectra to BC03 models and 
assume a Kroupa IMF. We only make comparisons for 
sources with PCA “warning” values set to zero. The 
PCA mass estimates suffer from aperture biases, which 
are increasingly important at lower redshifts, and from 
biases induced by low S/N CMASS spectra which be¬ 
come significant at higher redshifts. Unique to any of 
the estimates compared here, the Wisconsin masses in¬ 
clude priors on the presence of stochastic burst episodes 
(parameterized by the burst amplitude or mass fraction 
and a burst duration time) of constant star formation 
and the possibility of random, abrupt truncations in the 
star formation history. 

The Granada estimates are based on SED fits to the 
ES PS models (ESPS: Flexible Stellar Population Synthe¬ 
sis IConrov et al.l I2009L1201011 and a Bayesian approach 
similar to that adopted for our fiducial estimates. Four 
sets of priors are adopted which change the nature of 
ranges allowed for metallicity, star formation history (pa¬ 
rameterized with T models), age (less than the cosmic 
age), and dust content. The “Early” priors restrict to 
non-star-forming templates. “Wide” priors allow mod¬ 
els with more recent star formation. Both of these are 
available with and without potential dust extinction. We 
compare our results to Early (with dust). Early (no dust) 
which is similar to the M13 template, and Wide with dust 
(which is most similar to our fiducial priors). There are 
no burst models in the Granada templates. Both median 
and peak values of the M* posterior are provided. We 
use the median reported values, but choosing “mean” 
values does not make a significant difference, although 
the scatter increases for the “wide” set of priors. 

Gomparisons of the Portsmouth and PGA masses to 
our near-IR S82-MGC estimates are shown both as a func¬ 
tion of near-IR M* and redshift in Figure (TH Figure 
dn] presents the same comparisons against observed mag¬ 
nitude and the 61000 birth parameter (Section 18. 1|) . A 
similar set of comparisons against the Granada masses is 
displayed in Figures ITHl and ITTI 

8.2.3. M* comparison analysis 
https://www.sdss3.org/drlO/spectro/galaxy_portsmouth.php 
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Fig. 14.— Comparison of public BOSS M* estimates from [M^araston et al.l II2013I) and IChen et alj 1120121^ to those presented here (and 
labeled “near-IR”). Every panel plots the difference in logM* against near-IR M* in the left column, and against redshift in the right 
column. The Scr-clipped averages (light grey circles) and Itr standard deviation (gold lines) of the difference distributions are over-plotted. 
Shaded contours with levels separated by 0.3 dex increases in data density are displayed. The comparison is restricted to logM*/M 0 > 11.2 
for the right-hand panels. Overall normalization differences at the 0.1 dex level are expected. Of greater importance are possible systematic 
trends that may confuse evolutionary interpretations. Differences in M* estimates display little or no systematic trends with M*, but 
biases, likely resulting from a number of factors, are more apparent as a function of redshift. See text for discussion. 
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Fig. 15.— Mass comparisons as in Figure ^4] plotted now as a function of observed magnitude and the 6 iooo parameter, a proxy for the 
degree of recent s tar formation less sensitive to dust than an NUV-optical color. Some of the redshift-dependent trends in the M* offsets 
evident in Figure fT^ may result from a more fundamental dependence on magnitude. The comparisons reveal little dependence on 61000 ? 
despite the fact that the near-IR S82-MGC M* values are less sensitive to biases in MjL owing to recent star formation. The lack of any 
clear dependence on 61000 in the Portsmouth Passive comparison is particularly striking. Galaxies with higher 61000 values tend to be at 
greater redshifts, thus explaining the 61000 -dependent trend in the Wisconsin comparison. 
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Fig. 16.— Mass comparisons, as in Figure fT^ to the Granada M* estimates. We consider the median value of the Granada M* posterior 
distribution. The comparisons show relatively tight scatter but strong redshift-dependent trends that may be ind uced by a discrepancy in 
luminosity normalization, as evident in the more uniform trend with observed magnitude presented in Figure fTTI 
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CModelMag (i-band) Log B1000 

Fig. 17.— Granada M* comparisons as a function of observed magnitude and the 6iooo parameter for the Granada Wide estimates, 
which are the most similar in adopted priors to the S82-MGC measurements. The relative increase in the Granada M* at faint magnitudes 
continues smoothly across the full magnitude range, suggesting a difference in luminosity normalization with the S82-MGC estimates. 
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Despite different prior assumptions compared to the 
S82-MGC M* estimates and the use of different obser¬ 
vational data sets, Figures [14] through [T7| demonstrate 
agreement between the S82-MGC, Portsmouth, Wiscon¬ 
sin, and Granada estimates at the 0.1 dex level, with 
scatter in the comparisons of 0.2-0.3 dex that is consis¬ 
tent with the expected 0.1-0.2 dex uncertainties asso¬ 
ciated with any one estimator. As we explore possible 
trends below, it is important to emphasize that no single 
estimator can be considered “truth.” The S82-MGC esti¬ 
mates, however, benefit from more precise SEDs (thanks 
to the deeper Coadd photometry) and, because they are 
based on the near-IR, are less subject to potential biases 
from recent star formation. We emphasize that the com¬ 
parisons discussed here are restricted to bright, BOSS 
galaxies only. An important advantage of the S82-MGG 
is the ability to characterize galaxies fainter than BOSS 
and outside the BOSS selection boundaries in order to 
build M*-complete samples. 

In general, overall systematic offsets between M* es¬ 
timates are tolerable (and expected), but as shown in 
Paper III, the precision afforded by data sets such as S82- 
MGC now require understanding subtle dependencies in 
offsets at the 0.1 dex level. As a function of redshift or 
galaxy type, systematic biases or changes in the scatter 
at this level can strongly impact conclusions on evolu¬ 
tionary trends. 

Focusing first on the Portsmouth Passive masses, we 
see a relatively tight relation (upper panels in both fig¬ 
ures) against the S82-MGC estimates with a hint of in¬ 
creasing scatter at high 6ioqo values as might be expected 
if the passive template breaks down for galaxies with 
more recent star formation. The Portsmouth Best esti¬ 
mates feature more scatter and structure, perhaps reveal¬ 
ing how the combination of star-forming templates with 
a passive template produces a degree of disjoint structure 
in some regions of parameter space. 

The Wisconsin PCA masses feature a tight primary 
relation with the S82-MGC estimates plus a secondary 
distribution that scatters downward in the lower pan¬ 
els of Figure [m Figure [T5l Ibottom-left panel) suggests 
this scatter towards larger PCA M* values derives from 
i > 19 sources wh ere the BOSS spectral S/N drops and 
iChen et al.l (|2012h report biases in their own comparisons 
to SED-based M* estimates (although the sense of the 
bias appears to be reversed). A similar bias in the reverse 
direction was reported in M13. A gentle dependence on 
^1000 (roughly 0.1 dex across the sample) in the Wiscon¬ 
sin comparison is also visible, highlighting the effect of 
different priors in treating galaxies with more recent star 
formation (and higher dust extinction). Because the frac¬ 
tion of galaxies with higher feiooo values increases with 
redshift, this trend may underly the redshift dependence 
in Wisconsin M* offset at z > 0.4 (Figure [HI bottom- 
right panel). As the redshift decreases below z ^ 0.4, 
we speculate that aperture effects may drive the Wis¬ 
consin masses to lower values compared to the S82-MGC 
estimates. 

Turning now to the Granada masses. Figures ITSlandlTTl 
show some of the tightest relations (especially when plot¬ 
ted as a function of redshift, right-hand panels in Figure 
o compared to the S82-MGC estimates. This behavior 
may not be surprising given the similarity in priors used 
by Granada compared to the S82-MGC. The level of scat¬ 


ter is comparable to the z < 0.4 regime of the Wisconsin 
PGA comparison, which is encouraging because the Wis¬ 
consin group included bursts and truncated SFHs, but 
otherwise similar ranges for more continuous SFHs and 
stellar population parameters. 

Despite the tighter scatter, there are strong redshift- 
dependent trends in the Granada comparisons with z ~ 
0.6 galaxies about 0.3 dex more massive compared to 
the S82-MGC estimates than those at z ^ 0. This trend 
flattens somewhat at z > 0.4, but appears to continue 
when plotted as a function of observed magnitude (Fig¬ 
ure dZl left panel) suggesting a discrepancy in the over¬ 
all normalization of the luminosity. Understanding this 
discrepancy will be the subject of future work, but the 
lack of systematic trends as a function of M* or 61000 is 
nonetheless encouraging. 

9. SUMMARY AND CONCLUSIONS 

In this first paper in a series studying massive galaxies 
in Stripe 82, we have presented a new, publicly available 
compilation of Stripe 82 data sets, the Stripe 82 Mas¬ 
sive Galaxy Catalog (s82-mgc), designed to enable M*- 
limited studies of massive galaxy evolution since z ^ 0.7. 
The catalog includes 9-band ugrizYJHK photometry 
obtained by matching the SDSS Coadd and UKIDSS- 
LAS photometric catalogs using the SYNMAG package. 
Exploiting over 149,439 spectroscopic redshifts, we have 
assembled and tested a set of photometric redshifts from 
a variety of sources that are needed to complement spec- 
zs in order to build M, complete samples. With this red¬ 
shift information, we have derived a new set of near-IR 
based M* estimates and included them in the S82-MGC. 

A number of additional steps have been taken to make 
the S82-MGC ready for population studies. We have used 
the near-IR photometry to improve the star-galaxy sep¬ 
aration and built new estimates of UKIDSS total mag¬ 
nitudes that address the significant limitations due to 
blends in the public UKIDSS-LAS catalogs. We have also 
carefully defined the survey geometry and local depth 
of the combined data set, designing new custom masks 
in order to reject problematic regions in both the SDSS 
Coadd and UKIDSS imaging bands. 

Applying these rejection masks and a set of near-IR 
magnitude limits, we have constructed a M*-limited sub¬ 
sample of the S82-MGC called ukwide that spans 139.4 
deg^ and contains 41,770 galaxies with logM*/M0 > 
11.2 to z Ri 0.7, roughly 45% of which boasts a spectro¬ 
scopic redshift. This is the largest near-IR selected and 
M, complete sample of galaxies beyond z ^ 0.1 assem¬ 
bled to date. Paper II exploits this sample to study the 
completeness function in the BOSS samples. 

Finally, we have presented comparisons between the 
near-IR based M* estimates in the S82-MGC and previ¬ 
ous, publicly available estimates for BOSS galaxies. 
We analyze various systematic trends that are apparent, 
and find generally good agreement at the 0.1 dex level. 
However, as demonstrated in our forthcoming analysis of 
galaxy M, functions (Paper HI), systematics at this level 
are now a significant limitation in our ability to measure 
high precision evolution in the massive galaxy popula¬ 
tion. Encouragingly for the BOSS sample, the addition 
of near-IR data has only a mild impact on M* estimates, 
primarily affecting galaxies with more recent star forma¬ 
tion. 
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APPENDIX 
UKIDSS QUERY 


SElECr s . sourcelD, s. RA, s . dec , s.mergedClassStat, s.mergedClass, 
s . pStar , s . pGalaxy , s.pNoise, s.pSaturated, s. eBV, s . aY, s . aj , s . aH, s . aK, 

s . yPetroMag , s . yPetroMagErr , s . yHallMag , s . yHallMagErr , ydtcn . IsoMag AS YIsoMag , s . yppErrBits , 
s . yAperMagS , s . yAperMagSErr , s . yAperMag4 , s . yAperMag4Err , s . yAperMagd , s . yAperMagGErr , 
ymd. aperCorS AS yaperCorS , ymd. aperCor4 AS yaperCor4 , ymd. aperCorG AS yaperCorG , 

ymd. seeing AS yseeing , ymd. skyCorrCat as ySkyCorrCat , 

ymd. photZPCat as yPhotZPCat, ymd. photZPErrCat as yPhotZPErrCat, 

ydtcn . aperMagl AS yaperMagl , ydtcn . AperMaglErr AS yAperMaglErr , ymd.aperCorl as yaperCorl , 
ydtcn . aperMag2 AS yaperMag2 , ydtcn . AperMag2Err AS yAperMag2Err , ymd.aperCor2 as yaperCor2 , 
ydtcn . aperMagS AS yaperMagS , ydtcn . AperMagSErr AS yAperMagGErr , ymd.aperCorG as yaperCorS , 
ydtcn. ell AS yEll , ydtcn . pHeight AS ypHeight , ydtcn . pHeightErr AS ypHeightErr , 

s.j_lPetroMag, s.j.lPetroMagErr, s.j_lHallMag, s.j_lHallMagErr, j_ldtcn. IsoMag AS J_1 IsoMag , 
s . j.lppErrBits , 

s . j_lAperMag3 , s . j _1 AperMagSErr , s . j_lAperMag4 , s . j_lAperMag4Err , s . j.lAperMagG , s . j _1 AperMagGErr , 
j_lmd . aperCorS AS j_laperCor3 , j_lmd . aperCor4 AS j_laperCor4, j_lmd . aperCorG AS j.laperCorG , 

j_lmd. seeing AS j.lseeing , j_lmd . skyCorrCat as j_lSkyCorrCat , 

j_lmd . photZPCat as j_lPhotZPCat , j_lmd . photZPErrCat as j_lPhotZPErrCat , 

j.ldtcn . aperMagl AS j.laperMagl , j _ldtcn . AperMaglErr AS j.lAperMaglErr , j_lmd . aperCorl as j.laperCorl 

j.ldtcn . aperMag2 AS j_laperMag2 , j _ldtcn . AperMag2Err AS j_lAperMag2Err , j_lmd . aperCor2 as j_laperCor2 

j_ 1 dtcn . aperMagG AS j_laperMagG , j _ldtcn . AperMagGErr AS j_lAperMagGErr , j_lmd . aperCorG as j_laperCorG 

j_ldtcn. ell AS j_lEll , j_ldtcn . pHeight AS j_lpHeight , j_ldtcn. pHeight Err AS j_l pHeightErr , 


s . hPetroMag , s . hPetroMagErr , s . hHallMag , s . hHallMagErr , hdtcn . IsoMag AS HIsoMag , s . hppErrBits , 
s . hAperMag3 , s . hAperMag3Err , s . hAperMag4 , s . hAperMag4Err , s . hAperMagG , s . hAperMagGErr , 
hmd.aperCor3 AS haperCor3 , hmd.aperCor4 AS haperCor4 , hmd. aperCorG AS haperCorG , 

hmd. seeing AS hseeing , hmd. skyCorrCat as hSkyCorrCat , 

hmd. photZPCat as hPhotZPCat , hmd. photZPErrCat as hPhotZPErrCat, 

hdtcn . aperMagl AS haperMagl , hdtcn . AperMaglErr AS hAperMaglErr , hmd. aperCorl as haperCorl , 

hdtcn . aperMag2 AS haperMag2 , hdtcn . AperMag2Err AS hAperMag2Err , hmd.aperCor2 as haperCor2 , 

hdtcn . aperMagG AS haperMagG , hdtcn . AperMagGErr AS hAperMagGErr, hmd. aperCorG as haperCorG, 

hdtcn. ell AS hEll , hdtcn . pHeight AS hpHeight , hdtcn . pHeightErr AS hpHeightErr , 

s . kPetroMag , s . kPetroMagErr , s . kHallMag , s . kHallMagErr , kdtcn . IsoMag AS KIsoMag , s . kppErrBits , 
s . kAperMag3 , s . kAperMag3Err , s . kAperMag4 , s . kAperMag4Err , s . kAperMagG , s . kAperMagGErr , 
kmd.aperCor3 AS kaperCor3 , kmd.aperCor4 AS kaperCor4 , kmd. aperCorG AS kaperCorG , 


kmd. seeing AS kseeing , kmd. skyCorrCat as kSkyCorrCat , 

kmd. photZPCat as kPhotZPCat, kmd. photZPErrCat as kPhotZPErrCat, 

kdtcn . aperMagl AS kaperMagl , kdtcn . AperMaglErr AS kAperMaglErr , kmd. aperCorl as kaperCorl , 

kdtcn . aperMag2 AS kaperMag2 , kdtcn . AperMag2Err AS kAperMag2Err , kmd.aperCor2 as kaperCor2 , 

kdtcn . aperMagG AS kaperMagG , kdtcn . AperMagGErr AS kAperMagGErr, kmd. aperCorG as kaperCorG, 

kdtcn. ell AS kEll , kdtcn . pHeight AS kpHeight , kdtcn . pHeightErr AS kpHeightErr , 


1 . ymfID as ymfID , l.yeNum as yeNum, l.j.lmflD as j.lmflD , l.j.leNum as j.leNum , l.j_2mfID as j_2mfID , 
l.j_2eNum as j_2eNum, 1 . hmfID as hmfID, 1 . heNum as heNum, 1 . kmfID as kmfID , l.keNum as keNum 


ERCM 

lasSource AS s , lasMergeLog AS 1 , 

MultiframeDetector AS ymd, MultiframeDetector AS j_lmd , 
MultiframeDetector AS hmd, MultiframeDetector AS kmd, 


MultiframeDetector AS j_2md , 


lasDetection AS ydtcn 
lasDetection AS hdtcn 

WHERE 

s . frameSetID = 1 . frameSetID AND 
1 . ymfID = ymd. multiframelD AND 
1 .yeNum = ymd.extNum AND 
l.j_lmfID = j_lmd. multiframelD AND 
l.j.leNum = j_lmd . extNum AND 
l.j_2mfID = j_2md. multiframelD AND 
l.j_2eNum = j_2md . extNum AND 
1 . hmfID = hmd. multiframelD AND 
1 . heNum = hmd. extNum AND 
1 . kmfID = kmd. multiframelD AND 
1 .keNum = kmd. extNum AND 


lasDetection AS j.ldtcn , lasDetection AS j_2dtcn , 
lasDetection AS kdtcn 


1 . ymfID = ydtcn . multiframelD AND 
1 .yeNum = ydtcn .extNum AND 
s . ySeqNum = ydtcn . SeqNum AND 
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l.j.lmflD = j _1 dt cn . multiframelD AND 
l.j_leNum = j _1 dt cn . extNum AND 
s.j.lSeqNum = j _ldtcn . SeqNum AND 

l.j_2mfID = j _2dtcn . multiframelD AND 
l.j_2eNum = j _2dtcn . extNum AND 
s.j_2SeqNum = j _2dtcn . SeqNum AND 

1 . hmfID = hdtcn . multiframelD AND 
1 .heNum = hdtcn .extNum AND 
s . hSeqNum = hdtcn . SeqNum AND 

1 . kmfID = kdtcn . multiframelD AND 
1 .keNum = kdtcn .extNum AND 
s . kSeqNum = kdtcn . SeqNum AND 

(s.priOrSec <= 0 OR s . priOrSec = s . frameSetID) AND 
s.dec > —1.3 AND s . dec < 1.3 AND ( s . r a between 0 AND 30) 


S82-MGC DATA PRODUCTS 

Several types of s82-MGC data products are made publicly available at s82mgc.massivegalaxies.com. Documen¬ 
tation on the website provides the most up-to-date description of these products. We summarize what is available 
below. In some cases to make file sizes more manageable, we have divided Stripe 82 into east and west sections by 
splitting at Q!J2 ooo = 0 . 

Parent Catalogs 

Parent photometric catalogs in FITS table format for the east and west side of Stripe 82 are provided from queries of 
the SDSS Coadd and UKIDSS LAS databases. The Coadd parent catalogs have filenames beginning with “S82coadd_” 
and contain the Coadd photometry and SDSS profile fitting results. The UKIDSS catalogs have filenames beginning 
with “lasJDR8_” and contain the UKIDSS photometry, error flags, seeing, aperture corrections, and other information 
from the UKIDSS LAS database. 

The Stripe 82 Massive Galaxy Catalog 

The heart of the S82-MGC itself, with a filename starting with “pcatd_”, is selected on the SDSS Coadd parent 
catalog to include Coadd sources with measured r-band profile ht information (in practice sources with devAB_R and 
ExpABJl > 0.01 are selected). The peat table includes cross-referencing information for the UKIDSS catalog, matched 
SDSS-kUKIDSS SYNMAG photometry, and corrected UKIDSS total magnitudes, among other selected information 
from both parent catalogs. The peat FITS table is a single file that contains 15,342,585 with 92 tags for a total size 
of 5.5 CB. 

Galaxy Properties 

Additional FITS tables, matched one-to-one to the peat, provide estimates of galaxy properties for sources in the 
S82-MGC. These include spectroscopic redshifts from SDSS, VVDS, and DEEP2 and all photometric redshifts discussed 
in Section [H The S82-MGC fiducial M* estimates as well as KCORRECT absolute magnitudes and &iooo values are 
also provided. Please see the website for details. 

The UKWIDE sample 

For those interested in working with the M*-limited ukwide galaxy sample, smaller sub-catalogs of the peat and 
galaxy property tables are provided. Some software tools needed to construct ukwide, e.g., the star-galaxy separation, 
are also included. 

Survey Footprint 

The UKIDSS depth information and all mangle polygon files describing the geometric layout of the different survey 
components as well as the rejection masks are provided and described on the website. 
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